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Systems in which the dipolar energy dominates the magnetic interaction, and the crystal field generates strong 
anisotropy favoring the longitudinal interaction terms, are considered. Such systems in external magnetic field 
are expected to be a good experimental realization of the transverse field Ising model. With random interactions 
this model yields a spin glass to paramagnet phase transition as function of the transverse field. Here we 
show that the off-diagonal dipolar interaction, although effectively reduced, destroys the spin glass order at any 
finite transverse field. Moreover, the resulting correlation length is shown to be small near the crossover to 
the paramagnetic phase, in agreement with the behavior of the nonlinear susceptibility in the experiments on 
LiHoxYi_xF4. Thus, we argue that the in these experiments a cross-over to the paramagnetic phase, and not 
quantum criticality, was observed. 



The study of quantum phase transitions (QPT) is of prime 
recent interest, as it is beUeved that the understanding of the 
physics at the vicinity of quantum critical points will shed 
light on some of the most interesting problems in condensed 
matter physics, such as the metal insulator transition, super- 
conducting insulator transition and high temperature super- 
conductivity. Quantum magnets, and specifically their mod- 
eling by the transverse field Ising model (TFIM) 
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are a particularly good laboratory to study QPT, as this model 
is rich enough to capture the interesting physics of QPT, yet 
simple enough to allow theoretical treatment. Experimentally, 
much effort was invested to realize the TFIM, and maybe the 
best realization is in anisotropic dipolar systems. In these sys- 
tems the dipolar energy dominates the spin-spin interaction, 
and the crystal field generates strong anisotropy resulting in a 
ground state Ising like doublet for the single spins and an ef- 
fective reduction of all but the longitudinal interaction terms. 

Indeed, LiHoxYi_xF4 with x= 1 was shown jl|l to exhibit 
a ferromagnetic to paramagnetic (PM) transition as function 
of transverse field Ht and temperature T. As x is reduced, the 
randomness in the position of the magnetic Ho atoms results 
in frustratioii, and for x= 0.167 a spin-glass (SG) phase was 
observed |^|3|]. Furthermore, applying a transverse magnetic 
field induces quantum fluctuations, leading to a PM phase 
at large fields. Thus, this compound is considered to be the 
archetypal experimental realization of a quantum SG lAtSI- 

In this letter we show that for anisotropic dipolar glasses 
in general, and for the LiHoxYi_xF4 compound in particular, 
the off-diagonal terms of the dipolar interaction, albeit effec- 
tively reduced, qualitatively change the physics of the prob- 
lem. In particular, in the presence of a transverse field the 
off-diagonal dipolar terms reduce the symmetry of the system 
in comparison to the TFIM, and render the latter inadequate 
in studying the system. A proper treatment of the off-diagonal 
dipolar terms results in the absence of long-range SG order at 
any finite Ht, and a reduction of the cusp in the non-linear sus- 
ceptibility at the crossover to the PM phase. Thus, we argue 
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FIG. 1: Schematic picture for the T — behavior of a dipolar Ising 
SG in a transverse field Ht. The typical size ^ of a SG ordered do- 
main (depicted below the x-axis) decreases with Ht, with a critical 
exponent f (see text). At large enough Ht the system becomes PM, 
via a crossover and not a quantum critical point. 



that the experimental line drawn at the peak values of the non- 
linear susceptibilitylsl is not a phase transition line. Except 
for the point at Ht = 0, this line corresponds to a cross-over 
between a paramagnet to a phase we denote a "quasi spin- 
glass" (QSG). In this phase the system separates into domains 
within which the random ordering of the spins is maintained. 
These domains have a typical size i{Ht) which dictates the 
correlation length in the system, and its dependence on Ht is 
given by the critical exponent v calculated below. The domain 
structure is maintained until the crossover field, where fluc- 
tuations between the relevant Ising like states dominate and 
the system becomes PM (see Fig^for a schematic picture at 
T = 0). This crossover is expressed as a cusp in the non-linear 
susceptibility. Importantly, the reduction of ^ with increasing 
Ht results in the corresponding reduction of the cusp in the 
nonUnear susceptibility, explaining the peculiar experimental 
resultjsl where the cusp is reduced with decreasing T. Inter- 
estingly, we show below that at T = the crossover takes 
place at a value of Ht which corresponds to ^ w 1, and there- 
fore to a complete absence of a cusp of the non-linear suscep- 
tibility at r = 0, as can be inferred from the experiment^. 
Theoretical considerations. — Our analysis below is valid 
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both specifically to the LiHoxYi_xF4 system, as we further 
comment on below, as well as to any anisotropic dipolar sys- 
tem. The only requirement is that the single spins have a 
ground state Ising like doublet, with a large energy separa- 
tion to the excited states. To emphasize the generality of our 
approach we consider the following spin-s Hamiltonian 
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Here i,j denote the positions of the spins, randomly diluted 
on some lattice, denotes the dipolar interaction (a, /? = 
X, y, z), and > is the anisotropy constant due to the crys- 
tal field. For Ht = the GS of a single spin is doubly de- 
generate with §2 = ±s and zero energy. The corresponding 
states are denoted Its) and \ls). The first excited states have 
Sz = ±(s — 1) and energy ilo = (2s — l)D. Throughout the 
paper we assume that Oo ^ fJ-sHt, Vmax where Vmax is the 
largest dipolar energy between two spins in the system. We 
now define H = Ti.\\ + 'H±_ such that 

W,l = -DY,[{Stf - J E WStS^, (3) 



and 
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^^ = ^^E E vt,'sTs|~^JiBH,Y,st. (4) 

We assume that the dilution is such that 7i|| is equivalent to 
the classical Ising model with random interactions and ex- 
hibits a SG phase at low temperature. As this classical dipo- 
lar Ising SG is__equivalent to the short range Edwards An- 
(Eq.([0 with random nearest neighbor 
]), our analysis is done within the scaling 
("droplet") picture l9|] which accounts for its behavior at large 
sizes. The GS of 7i || is then two fold degenerate 0| with states 
I V'o), I i'o), which are related by Sz ^ —Sz symmetry, and 
in which each spin is in either state \]s) or \[s)- Importantly, 
adding a transverse field term preserves the above symmetry, 
and therefore the TFIM is the archetypal model for the quan- 
tum SG phase. However, when adding Tij^ which includes the 
off-diagonal dipolar terms, this symmetry is not preserved. 
The GS degeneracy breaks, and the system gains energy by 
choosing locally a state similar to | -0o) or | ijjo) according to 
which optimizes the energy gain due to Tii^. 

Following the scaling picture of Fisher and Huse |@| and 
using an Imry-Ma lHoll like argument we calculate this energy 
gain, i.e. the energy to flip a droplet of size L having N ^ L'^ 
spins, due to the addition of H±. This energy gain (see Eq.(|9} 
below) is then compared with the energy cost due to the do- 
main wall formation, and the correlation length is obtained 
(II U . Although, as is shown below, one can define an effective 
longitudinal random field at each site and use a direct anal- 
ogy to the Ising SG in a random field, we would proceed by 
calculating directly the energy gain. 
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The addition of 'Hj_ to TC\\ changes | ipo) ~^ | "0) ™d | V'o) — > 
I with energies and E^/ respectively. The energy the 
system gains by choosing locally the lowest energy state is 
6E = \E^ — E,i,i I which we now calculate. In second order 



(^2) perturbation theory E^, ~ E^ 



E^p'' where 
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Here we used the fact that the only relevant excited states are 
those in which one spin changes its state from Sz ~ ±s to 
Sz = ±(s — !)■ Therefore the energy of all relevant excited 
states is flo in leading order, and the sum over the excited 
states can be taken out as the identity operator. A similar 
equation holds for E^' . One can show that the terms with 
even powers of Ht are equal for E^ and E^/ , while the term 
linear in Ht is equal in magnitude but has opposite signs for 
Eji, and E^i 0. Using the fact that since k ^ / the operators 
commute we obtain 
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and therefore 
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where we define hf = J2k ^ki i'^k) effective transverse 
magnetic field at site i. For each i all the V^i's are small except 
the few for which the sites i and k are spatially close. Due to 
the randomness of the sign, retaining for each i the term with 
the largest absolute value, denoted Vi, gives a good approxi- 
mation for 6E up to a numerical factor c of order unity. Since 
Vi is random in sign, the average energy gained by flipping a 
droplet of N spins is given by 



{6E) = 
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and we choose 



where V is the average magnitude of 

I -0o) and I i'o) such that 6E > 0. 

The above result is central to our analysis, and in order 
to check our approximation of randomness leading to it we 
calculated the gap between the GS and the first excited state 
numerically using Lanczos exact diagonalization (ED) fllll . 
We consider system sizes in the regime where they are much 
smaller than ^. This is important for our calculation, since 
then we are dealing with single domains, and therefore the 
two lowest states correspond to | ■0) , | V'')' and the gap to 6E. 
To reproduce the experimental situation, we focus on three di- 
mensional finite size clusters, randomly distributing N spins 
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at the rare earth sites of the LiHoxYi_xF4 lattice. Since we 
are interested in small fields, it is sufficient to use s ~ 1 parti- 
cles with on-site anisotropy JIq. We therefore study the spin-1 
version of H = Hn + H , : 



(10) 



with 6 w 0.2 in 3 dimensions fl4l[l5ll . Furthermore, under 
quite general conditions Fisher and Huse argued 01 that 9 < 
(d— 1)/2. For the dipolar Ising SG we expect the same scaling 
behavior with a similar exponent 9^^^ 6 to hold 0,0,1^, and 
the energy of flipping a domain is therefore « Vs^L^''. As a 
result, for L such that (s'^ ^BHt VVN)/no > Vs^L'^" it will 
be preferable for domains to choose their state between | ip) 
and I ip') as the one that locally minimizes £'^^ . This results 
in a finite correlation length 



Here and below all energy scales are expressed in units of 
the typical n.n. dipolar energy Vq. We fixed the dilution to 
a constant x = 3/16 = 18.75% by using 2 x 2 x N/3 unit 
cells, N being the total number of spins (there are 4 rare earth 
sites per unit cell). For this dilution we find that V = 0.8. 
Thus the Lanczos ED have been performed in the full d = 
[2s + 1]^ dimensional Hilbert space for = 6, 9 and 12 
s = 1 spins fTiil . Then, the gaps have been computed for 
10, 000 independent random samples for each size. 
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FIG. 2: Distribution P(^ln 

ED data obtained for the spin-1 Hamiltonian <10> with fio = 50 
and usHt — 0.5, have been collected over 10,000 random samples 
with the LiHoxYi_xF4 structure and a dilution x — 18.75%. Three 
different sizes have been used: A'^ = 6, 9 and 12, as indicated on the 
plot. Inset: Linear dependence of the disorder average gap (SE) vs 
\/iV. The dashed line is a fit of the form a\fN with a ~ 0.008. 



In Fig. 12] we present the numerical results obtained in the 
perturbative regime (i.e. D,o ^ /^B-^t and {ii^Ht)^ /VLo <C 
Vraax)- The ^/N Scaling of 5E as stated in Eq. (|9} is cjearly 



demonstrated, as we found a very good data collapse yJI for 
the distribution of In The inset of Fig.|2]also shows that 

the disorder average gap {5E) = a^fN. Confronting the nu- 
merical estimate obtained for the prefactor a with Eq. (|9j, 
we get c ~ 1. We have also checked the scaling of the gap 
with Ht/fto, where for several combinations of {ilo,Ht), we 
obtain an excellent collapse of the data by rescaling 5E 

Ht/na Lai- 

In order to obtain the correlation length of the system, i.e. 
the typical domain size, we have to compare the domain's en- 
ergy gain |9] to the energy cost due to the formation of a do- 
main wall. For the short range Ising SG this energy is cx 
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For the Ising SG in longitudinal field it was argued 0, [l^ 
and then verified experimentally 

US El El and numeri- 
cally i2(i[Elll that there is no de Almeida Thouless line IH, 
and no SG phase at any finite field. At Ht <C fio/MB our 
system is equivalent to the above model in small longitudinal 
fields, and we thus argue that there is no SG phase at any finite 
transverse field when the interaction is dipolar, and as Ht ^ 
the correlation length diverges with the same form H of the 
critical exponent ly = f^^/2)-0d ■ 

In our treatment the only dipolar terms we considered are 
the longitudinal and the zx terms. However, one can show that 
all the neglected terms [see the terms present in 7i_L, Eq. 
but not in Tij^, Eq. (BJl do not contribute to 5E in second order 
perturbation theory j7|. 

Interestingly, the two effects of the transverse magnetic 
field, i.e. inducing the crossover to the paramagnetic phase, 
and the reduction in ^ calculated above, behave very differ- 
ently as function of Ht- The former is dictated by fluctuations 
between the two single spin Ising ground states, which de- 
pend on Ht to a high power, of order s, and are practically 
negligible as long as Ht ^ 57o/mb- However, the fluctua- 
tions that dictate the reduction of ^ at low transverse flelds 
are between each single spin ground state and its flrst excited 
state at energy VIq. The latter depend on Ht to second order 
and result in a reduction of ^ which depends on 1/Ht to a 
power V close to unity. Therefore, the disordering of the SG 
order by Ht occurs in two stages. At low field domains of 
size ^ are formed, within which the GS is very similar to ei- 
ther of the two zero field SG ground states. At Ht ~ ilo/ Mb 
a crossover occurs where the order within each domain is de- 
stroyed, and each spin in the system fluctuates independently. 
Importantly, when reaching the crossover region at very low T 
one is already in the regime where ^ w 1 in units of inter-spin 
spacing, resulting in small features in the relevant suscepti- 
bilities, in agreement with experiment QS- We would like 
to emphasize that the understanding of the scenario above re- 
quires a model in which the large spins are considered and 
the anisotropy in explicitly taken into account. Indeed, the 
anisotropy energy JIq enters explicitly into Eas.(l9ll.( fTll . A 
presumably simpler model in which one treats spin-half par- 
ticles and models the effective reduction of the off-diagonal 
terms in the dipolar interaction by a multiplicative reduction 
factor will not be sufficient, since in such a model both the re- 
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duction in the correlation length and the crossover to the PM 
phase are induced by fluctuations between the Ising ground 
states, and therefore have the same scale in magnetic field. 

In addition to changing the symmetry of the system at 
Ht 7^ 0, resulting in the destruction of the SG phase and the 
QPT to the PM phase, the off-diagonal terms of the dipolar 
interaction also enhance the effective transverse field lEsll . Al- 
though, in principle, [see Eq.(|8j]is a random quantity, do- 
mains of size ^ choose to be in a state equivalent to | ip) or 
I ip') by the maximization of hf. As a result a net mag- 
netic field in the x direction, (hf) = {5E)/N oc is 
added to the external transverse field. As the crossover re- 
gion is approached ^ is small and the effective transverse field 
due to the off-diagonal dipolar interaction is significant. We 
thus give a precise physical origin to the conjecture made in 
Ref.Hl 

Our analysis above could equally be done by defining 5E in 
Eq. © as Y.k KiSl), where hi ^ {2s^Ji^Ht/n,) V^^. 
Using this definition one can make the analogy between the 
current problem to the Ising SG in random longitudinal field, 
as an alternative to the direct calculation of 5E performed 
above. 

Experimental consequences. — The crystal field Hamil- 
tonian in LiHoxYi_xF4 is different from the one given 
in Eq. (|2j. Furthermore, the strong hyperfine interactions 
strongly re-normalize the parameters of the TFIM, invalidat- 
ing the simple model in the electronic degrees of freedom lEsIl 
Still, for our purpose here an equivalent physical picture 
emerges; the two relevant (electro-nuclear) Ising states of 
each Ho ion couple very weakly at small transverse field, and 
the relevant excited states are at w lOK above the ground 
states Thus, the requirements for the validity of our theory 
given before Eq. (|2} are fulfilled. Our analysis and results 
[and in particular Eq. (IllH are therefore directly applicable to 
the SG experiments in the LiHoxYi_xF4 system 121 13|, with 
51o ~ lOK and suggest that LiHoxYi_xF4 is not a SG at any 
Ht ^ 0. Furthermore, the peculiar experimental result |y!] 
where the cusp in the nonlinear susceptibility is reduced with 
decreasing T is naturally explained; as T is reduced the 
crossover to the PM phase occurs at larger transverse fields. 
This results in smaller correlation length ^, and therefore a 
diminishing of the cusp in the nonlinear susceptibili ty ibj . In 
addition, the re-normalization of the effective spin 112 3ll spe- 
cific to the LiHoxYi_xF4 compound further reduces the non- 
linear susceptibility near the crossover 

From the experimental point of view our analysis changes 
the status of the field. The only claim for the observation of 
the QPT between the SG and PM phases was made in RefQ. 
As it is clear by our analysis here that it is a crossover and 
not a phase transition that was observed at low temperatures 
in the above experiment, an experimental observation of this 
QPT is still awaiting. Our analysis also points to the direction 
one should take in seeking such a QPT; systems in which SG 
order and quantum fluctuations compete, and either or both 
are controlable by a parameter which does not change the 
symmetry responsible for the GS degeneracy of the ordered 



state. An example would be the change, with applied pres- 
sure, of crystal field terms which induce quantum fluctuations 
between the Ising like doublet (such as (5^+5^) terms added 
to the Hamiltonian Q for integer spin systems). 

Recently there is an increasing experimental 
and numerical lE(i. ETI support for the validity of the droplet 
picture in describing short range Ising SG in general, and to 
its prediction 116(1 of the non-existence of a de Almeida- 
Thouless t22il in particular For the anisotropic dipolar sys- 
tems discussed here the crossover to the PM phase at Ht « 
^o/fJ-B is a result of quantum fluctuations, and there is no 
analog to the de Almeida-Thouless line. However, at Ht <C 
^o/fJ-B the system is equivalent to a classical Ising SG in a 
small random longitudinal field. Thus, the above numerical 
and experimental results 20, 2 IJ support the va- 

lidity of the droplet picture for dipolar Ising systems in small 
transverse field as well. Still, we believe that experiments that 
would directly observe whether dipolar Ising glasses in gen- 
eral and LiHoxYi_xF4 in particular have a SG phase at a fi- 
nite transverse magnetic field are of much interest, both as a 
verification of our results, and as an additional support for the 
droplet picture in general. 

Finally, our analysis is applicable also to any Ising SG 
where the dipolar interactions are present, even if the inter- 
action that governs the ordering is different. In such case the 

correlation length will be given by 0] w ^ JgHtV ) ' ^ ' ' 
where J is the strength of the dominant interaction. The qual- 
itative picture will remain similar, only now the size of the 

domains at the quantum crossover to the PM phase would be 

_ ^ j_^i/{3/2-e) 
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